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CN ■ Abstract 
> ■ 

■ We studied the variation of Newton's gravitational constant G, using 

' the experimental data and the scalar-tensor theory for gravitation. We 

, search for a fine tunning to match the datas and find possible scenarios 

^SJ ' for the history of this constant. 

o : 

m . 

O ■ 1 Introduction 
o '. 

Since the thirties, many cosmologists have been suggested theories that predict 

^ I the variation of Newton's gravitational " constant" . Among the more famous 

bj). are Dirac's Large Number Hypothesis [1], the Brans-Dicke theory [2] and more 

^ ' recently, some model of scalar-tensor theories (see the work by Barrow [3] . for 

. a review). The last one can be seen as a generalization of Brans-Dicke theory. 

' On the experimental hand, some measurements are performed to study this 

^ . variation; for instance, the spacecraft radar ranging [4] and lunar laser ranging 

that is: 



[5] . Both experiments give an upper limit for the ratio G/G 

G 



On despite of so especulative ideas, I have seen few numerical analisys 
about the evolution of the value Newton's gravitational constant. A work by 
Bronikov, Melnikov and Novello [6] is in this sense. They describe the upper 

limit to G/G . Barrow and Parsons [7] studied the behavior of cosmological 
models with varying G .They describe some analytical solutions to the field 
equations, using the conformal time coordinates and showed the dependence of 
the coupling constant lo with 4>/ 4>q in some particular cases of the models. 
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In the present work, we made a numerical analisys, based in a Scalar-Tensor 
theory, with variable w, studied by Barrow and others [3]. These numerical anal- 
isys, togheter with experimental data, could indicate the more correct theory, 
the initial conditions of the Universe and the evolution. The paper is divided in 
the following sequence: in next section, we review the model and isolated some 
variables that will be important in the sequence. In section 3, we show the nu- 
merical methods and hypotheses. In section 4 we show some graphics, the main 
point of this work. The last section is devoted to discussion and concluion. 



2 Relativistic scalar-tensor theories 

Newtonian gravitation permits us to ' write in' an explicit time variation of G 
without the need to satisfy any further constraint. However, in general relativity 
the geometrical structure of space-time is determined by the sources of mass- 
energy it contains and so there are further constraints to be satisfied. Suppose 
that we take Einstein's equation in the form (c= 1) 

Gl = 8nGT^ 

where and are the Einstein and energy-momentum tensors, as usual, but 
imagine that G = G{t). If we take a covariant divergence -a of this equation 
the left-hand vanishes because of the Bianchi identities, Tj^.^ = if energy- 
momentum conservation is assumed to hold, hence dG/dx'^ = always. In 
order to introduce a time- variation of G we need to derive the space and time 
variations from some scalar field, ijj which contributes a stress tensor T^{'4>) to 
the right=hand side of the gravitational field equations. 

We can express this structure by a choice of lagrangian, linear in the cur- 
vatura scalar R, that generaizes the Einstein-Hilbert lagrangian of general rel- 
ativity with, 

L = -/(V)i?+^a„Va> + 167rL„ (1) 

where is the lagrangian of matter fields. The choice of function f{tp) defines 
the theory. When i/j is constant this reduces, after rcscaling of coordinates, to 
the Einstein-Hilbert lagrangian of general relativity. The formulation introduced 
by Brans-Dicke [2] can be obtained from 1 by a non-linear tranformation of 
and /. Define a new scalar field ^ , and a new coupling constant Lo{(j)) by 



then 1 becames 



(2) 



-<j)R + + lewLm- (3) 



The Brans-Dicke theory arises as the special case 
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Brans — Dicke : t^((/>) = constant; f{ijj) oc ijj'^ . 



(4) 



The field equations that arise by varying the action associated with L in 3 
with respect to the metric, gab and (j) separately, gives the field equations, 



<l^,a<P,b - 1^9ab(l),i<P,i 



- <P ^ [<P,a;b - 9abO(j)] (5) 



[3 + 2a;(</.)]n</. = 87rT«-a;'(</))<A,, 



(6) 



where the energy-momentum tensor of the matter sources obeys the conserva- 
tion equation 



nab 



0. (7) 

These field equations reduce to those of general relativity when (and hence 
ui{(j)) ) is constant, in which case the Newtonian gravitational constant is defined 
hy G = cl>-\ 



2.1 The Basic hypothesis of our anaUsys 

In this work work, wc will not suppose any trial function io{(j)), hence wc will 
not propose any analitycal solution to the field equation; instead we calculated 
G{t), numerically, based in initial values of ^,<^ and of course uj{4>) (given by 
the field equations). We suppose that the scalar tensor theory can describe the 
present Universe quite well, so we have to consider the recent experimental data 
about it. We also assume that these data are hold, at least, during the last one 
billion of years. 

Clearly the equation of state of our work must be: 



p = 0, (8) 

that is a basic fact in the cosmology. 

Now is official, the Universe is flat. With this headlines, the Physicsweb 
site call the attention to work by Bernardis et al., in the Nature [8] , about the 
measures of the density of the Universe. This give us another information to 
work, the tricurvature: 



A; = (9) 

Last, but not least, is the observation that the Universe is expanding [9], not 
only expanding, but also accelerating, that is, the deceleration parameter q is: 

^ = -^ = -1 (10) 

where a is the scale factor, that implies : 
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a 

where H is the Huble parameter, and hence 

a = e"\ 

The field equations we have to work now are as follow: 
For the Robertson- Walker's metric, 

rfs^ = dt^ - a{tf [dx'^ + dy^ + dz^] 
we find, then {Lm = 0) : 

-.2 

H 



Sep 0^ ^ 6 



I + 3H(f) ■■ 



3 + 2a; 3 + 2a; 



p+3-(p + p)=0. 



One can isolate a;. 



u> = Q- 



-am 



3(j) 



and o), 



a; = 12. 



6^ 



2 STTPoe-^^* 



-H 



2 1 

3? 



{247rpo£ 



-3Ht 



(11) 

(12) 

(13) 

(14) 

(15) 
(16) 

(17) 



(18) 



Here Pq is an initial value for the density. Putting uj and uj in 15, we can 
isolate 0. This is an important point of this work, we use the diference finite 
method to calculated the field variations. We can check that: 



2H + Snpoe-^"^ - 7H(t> - 12(070) - 87rpoe=^^V(30) + 0iI/0 



and can use : 



1 + 2iI0/0 



^{t + At) = 0(i) + 0jAf 



(19) 
(20) 



In the same way: 



(t){t + At) = (f>(t) + ^At (21) 

and we have recorrence. We still need some information, the intial values for </> 
and 6 and will discuss it in the next section. 



3 The initial values problem. 

As we can see in 19, depends on the density pg, the Hubble parameter H, 
that are constants, (p depends also of ^ and (j) and the choice of these two 
intial values determines the evolution of Newton gravitational constant. We 
remember that: 

Hence u{^{t)) and ^{t) must result, when we arrive at the present Universe: 

G = 6.67 X lO-^dyn.cm^g-2 (23) 

Note that we start in the past and go to the present time, so we can have 
values of G in the past. 23 is a final value condition. The theory has two free 
parmeters, the initial values for (j) and (j), and could have also two constraint, G 
and G, but there are few measures of G, so still there are some freedom in the 
choice of initial values. At the end of this work, we show a result that consider 
the limit value for G. In the Brans-Dicke [2] theory, we have 

<^oci. (24) 

We use this fact as a guide to the choice of (p. For (j) we use the relation: 

(/) ~ 10"^. (25) 

n is an integer. If we write: 



(26) 



where 6 is a constant, each n value will determine the value of 6, (or we can also 
work in the opposite sense) in order to satisfy relation 23. This is the important 
point of this work, we have to choice the initial num,erical values for (j) and for (j). 
These values held in the past, one billion years ago. So the ratios 4>/4> discussed 
here must be considered in the past. 

Finally, for H we use the same value of [6] that is : 

H = 0.7x 10-^°yr-^ (27) 
Pq comes from [11] and it is : 
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(28) 



Of course, when considering the past, we have to correct this density. 

With these considerations, we have performed some possible situations, that 
are presented in the next section. 

4 Results. 

Wc studied some possible cases, using the datas and methods explained above, 
taking care that each step in our program means 10 millions of years. So in the 
figures below, 100 means one billion of years ago. 

The first case is the following: we have positive values for cj), in a way that: 

4> = 10-"?!>. (29) 

Wc performed our analisys for the case n = 2,3,4,6,8. In the figure 1 we 
present our results for the case n = 1, 2, 3. In our picture, the present time is 
on 0.0 and the past goes to right. 

In figure 2, we studied another values. We use n = 4, 6 and 8 and (/) is 
positive. In the figure 2 we can these results. 

As we can see, it's interesting to work the case n=3. So we studied the 
cases: n = 2.5, and the more interesting: 

^ = 2 X 10-^(t>. (30) 

In this case, the variation of G{t) is slow and the calculated value for u){(p) 
is bigger than 500. This is in agreement with obscrvacional data [10]. More 
exactly, a;((/>) w 889 (for the present time). This means that this is a good 
starting point to search the fine tuning of the theory with the experimental 
data. In the figure 3 wc have the results. 

For the sake of completeness, we also show the cases where we have a negative 

^. 

4.1 The fine tuning 

As a final result, after a search for the parameters that hold theory and expri- 
mental datas, we find that using: 

(j) « 0.92755/G 

as the intial value (one billion years ago), and (also for one billion years ago) 

</) = 1.754236908125 x Ur^cj) 
we obtain the one posible value for G, one billion of years ago: 
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Gpast ~ 7.54 X 

and the value for uj, in the past, was, according to our computation 

l^past ~ 8.87 

and the present value for w is, according to our program: 

w » 8.24 X 10^^ 

and finally, the present ratio: 

Q 

— « -5.77 X 

Gr 

Note that, because we have so small experimental value for this ratio, we have 
to use several digits in (j) . 

These are the initial results of our proposal of perform an numerical analisys 
of the consequences of the scalar tensor theories of gravitation. We hope that 
these informations can be used in another investigations of the cosmological 
parameters. 



5 Discussion and Conclusion. 

In this work, we consider the scalar-tensor theory of gravitation, which the 
has the proposal of study of the Neton's gravitational constant. We didn't 
sugcst any analitycal solution to the field equations. What wc have done is, 
based in some justified hypothesis, to calculated numerically the solutions for 
the fields {(j)), the coupling constant (w), and the main point, the values of 
G{t). On despite that the field equations are non-linear in (j) and (p wc have two 

constraints (G and G/G ) that filter the possibilities. It's useful to study if 

we can aply the uniquiness of solution under these conditions. 

We consider that numerical analisys presented is quite interesting, because 
shows the possible evolutions for the value of the Newton's gravitational con- 
stant G in a recent period of the Universe, based in a Scalar Tensor Theory of 
gravitation. The same idea can be used to make some predictions about future 
values of G. We believe that, if one consider our basic hypothesis valid, (fc = 0, 
accelerating Universe, p = Q) it's worthwhile to work to find the intial values 
of (j) and 4> that results in the present values of G and G (if the last one is 
measured). The numerical approach used here can be improved to work with 
a more accurate step in the program (less than 10 millions of years). Another 
point is that our computational program shows that is possible a change in the 
signal of G, so we may especulate that G can start to grow. 

We can also consider that, if the Newton's gravitational constant have 
changed in the last one billion of years, there are consequences and evidences 
that can be search in the geological or astronomical fields of knowledge, for 
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a comprccnsivc discussion about these implications, sec [11]. Eventually, the 
study of gravitational waves can give some information about the value of G in 
the past, or , in the opposite sense, these study can use some information about 
the numerical evolution of G, gravitational waves in the scalar-tensor theories 
are studied in [12]. 

As a final consideration, we remember that we use Lm = 0. It's possible 

propose many different forms for Lm- In our humble opinion, if we think that 
even the theorethical part is not solved yet, hence the numerical one, in a way 
that it 's possible to perform several approachs to the problem. 
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Captions 



Fig 1:G as a time function, when we consider the some ratios 
between (f) and 0, in this case 4'/4' = 10~", where n — l(solid line),2(long 
dashed line) or 3(shor dashed line) . It 's important to observe that 
in this picture, the time goes from the right to left, being 0.00 the 
present time and 100 one billion of years ago. The same notation is 
used in the another grafics. 



Fig 2: The case n=4, (solid) 6 (dashed) and 8. ^ is positive. We 
can see that n=6 and n=8producesalmost the same results, so one 
line superposes the another. We can also see that we have a growing 
G(t), and also note this is an acceleratting growing. Due this we will 
eliminate these results of a "fine tunning" analisys. 

Fig 3:A study for the cases with n=2.5(long dashed) and n=3, 

((/) = 2 X 10~^(/)) (short dashed). For (f) = lO^'^ct) (solid line) we can observe 
a minimum. We can see there are an interesting result and worked it 
as a starting point of a search for a better fitting of the theorethical 
peu-ameters. 

Fig 4:G(t) in the cases we have negative values for the intial value 
of 4). We use n=4, (solid line) , 6 (long dashed) and 8 (short dashed). 
Again one result superposes another. The cases n=3 and n=2 have 
resulted inconsistent values and oscilations. 



9 




Fig 1 
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